Dynamical Mean-Field Study of the Ferromagnetic Transition Temperature of a 
Two-Band Model for Colossal Magnetoresistance Materials 
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The ferromagnetic (FM) transition temperature (To) of a two-band Double- Exchange (DE) 
model for colossal magnetoresistance (CMR) materials is studied using dynamical mean-field theory 
(DMFT), in wide ranges of coupling constants, hopping parameters, and carrier densities. The 
results are shown to be in excellent agreement with Monte Carlo simulations. When the bands 
overlap, the value of Tc is found to be much larger than in the one-band case, for all values of the 
chemical potential within the energy overlap interval. A nonzero interband hopping produces an 
additional substantial increase of Tc, showing the importance of these nondiagonal terms, and the 
concomitant use of multiband models, to boost up the critical temperatures in DE-based theories. 
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PACS numbers: 75.47.Lx, 71.27. +a, 71.30.+h 

The study of the CMR rare-earth perovskites has re- 
cently attracted considerable attention due to the rich 
magnetic and structural transitions they display Q . The 
CMR effect, namely an extremely large drop in resis- 
tivity caused by a magnetic field, occurs at the tran- 
sition between a low-temperature FM-mctallic ground 
state and a high-temperature paramagnetic-insulating 
phase, i.e., near Tq. A basic aspect of CMR physics 
is that the mobile carriers (Mn e g -symmetric electrons) 
are strongly coupled ferromagnetically to localized spins 
(Mn ^-symmetric electrons). The core spin alignment 
influences on the electron motion leading to ferromag- 
netism. The basic model used to describe CMR mate- 
rials is the DE model Q, formulated with one localized 
spin per site coupled to mobile carriers via a Hund's cou- 
pling J much larger than the hopping amplitudes. While 
DE ideas explain qualitatively the existence of ferromag- 
netism, much of the CMR physics is still under debate 
0, • This is in part caused by the absence of fully reli- 
able many-body techniques to study the complicated DE 
models needed to describe these compounds. Only MC 
and static mean-field approximations have been applied 
to the realistic multiband problem the DMFT treat- 
ment being notoriously absent. Despite a considerable 
effort carried out within DMFT in the context of one- 
band models [f|, the realistic case of two active bands 
has received much less attention since its analysis is far 
more difficult. 

In this paper, the DMFT method is applied for the 
first time to the analysis of Tc of a two-band model for 
CMR materials. Its results are contrasted against MC 
simulations and found to be in good agreement. Besides 
the relevance of this combined DMFT+MC contribution, 
which is highly nontrivial technically as shown below, we 
also unveil here the key relevance of the interband hop- 
ping amplitudes to boost the value of the critical FM 
temperature. This effect was also not discussed in pre- 
vious literature, and may lead to creative procedures to 



further increase Tc in real materials. The DMFT ap- 
proach is general, with two s=\/2 active bands, arbitrary 
couplings, hoppings, and carrier densities p. 

The two-band model is based on the DE Hamiltonian: 
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where I, I (= 1, 2) are the e g -band indexes, i,j label the 
sites, ci t i a destroys an electron at site i in the band /, 
&l,i,pa = c\ i p((Tf3 a /2)ci ! i ta is the spin- operator of the 
mobile carrier at band / (<r = Pauli vector), a and (3 are 
spin indexes, Ji is the coupling between the core spin and 
the conduction electrons of band I, and Si is the localized 
spin of the magnetic ion at site i. For I = I we refer 
to tu as the band hopping (= t{). For I ^ I , the t n ' is 
referred to as the interband hopping (t u > =t l > l ). The two 
active bands couple through the simultaneous scattering 
of carriers with the same core spin, as well as through 
the exchange of carriers via the nondiagonal hopping. 
While the first coupling clearly causes an increase in Tc 
when the bands overlap within the same energy interval, 
namely when Ji are closed, the electron exchange among 
bands is also shown to induce higher Tc even when the 
Jz's are very different. 

Note that the addition of other important terms, such 
as the antiferromagnetic exchange Jaf between core 
spins and/or the cooperative Jahn- Teller phonons, will 
need a sophisticated "cluster" DMFT where at least some 
short-distance effects are considered. For this reason, this 
first study using the DMFT method on a two-band model 
focuses on the simplest case where only the FM phase is 
relevant. The study of the competition with other phases 
is left for future investigations. The model is studied be- 
low using DMFT and MC techniques. 

DMFT results: Within DMFT, the self-energy is lo- 
cal; S(p, iw„)— >£(iw n ) (uj n =(2n + l)irT are the Mat- 
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subara frequencies). Since £ is momentum indepen- 
dent, the information about the hopping of carriers 
on and off lattice sites is carried by the bare Green's 
function Q (iui n ). With the effective action 5 e //(m) 
defined by Qq quadratic in the Grassman variables, 
the full Green's function Q can be solved by integra- 
tion: (G(iu} n )) = ([Go (iu n ) + JSma] -1 ). The aver- 
age (X(m)) = J dfl m X(m)P(m) is over the orientations 
m of the local moment with the probability P(m). In 
the FM phase near T c , P(m) oc exp (— 3/3Mm), where 
/3 = 1/T, and M = (m z ) m is the local-order param- 
eter, nonzero only below Tq. If two bands are ac- 
tive, then one rewrites the average above for each band; 
(Ql(iuj n )) = ([Q^l(iu} n )+ JfSjma-]^ 1 ), and solve the equa- 
tions for the coupled Green's functions; (Gqj (iu} n )) = 

z n - tf (Gi(iuJn)} ~ t 2 u > (Gi'(iu n )), (I ^ on a Bethe lat- 
tice with a semicircular noninteracting Dcnsity-Of-States 
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(p = chemical potential). t U ' carries in Go,i the infor- 
mation about the second band I through Gi 1 . To find 
T c , we parametrize GqJ- Goj(iu n ) = [z n + Ri(iui n )]i + 
Qi(iu) n )a z , and linearize GqI with M. Up to first order: 
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where Bi(iuj n ) = z n + Ri(iuj n ). A similar equation for 
can be obtained by interchanging the band indexes I — ► I . 
Despite the extreme complexity of Eqs. (J5J), we managed 
to perform analytical calculations to decouple the equa- 
tions set and to obtain separate equations for Ri and R^ . 
The polynomial coefficients of the 9-th order equations 
for Ri and R^ which were solved numerically, are combi- 
nations of all parameters of the model, i.e., couplings and 
hoppings (not reproduced here because of their size). For 
t u > = 0, the 9-th order equations reduce to the 3-th order 
ones: Rf + 2z n R? + (z 2 n + tf- jf)Ri + z n tf = 0. At fj, = 
and with the substitution iuj n — > u>, Eq. (J2J gives the in- 
teracting electronic DOS; at T = 0. Due to having a 
nonzero t,f , besides its own parameters U and J;, DOS; 
depends also on the parameters characterizing the other 
band I , i.e., and Jf . This interplay of parameters de- 
scribing the coupling of bands leads to an increase in Tc 
even when the Jj's are very much different, i.e. when the 
bands occupy different energy intervals. The obtained 
equation for Qi(iu n ): 



Qi = % 



JiM 



Bf 



J? «' B 2 , - J?, 



+ tl 3{Bf 



2J 2 Qi 



+ t 2 , 



V V 



3{Bf, 



(3) 



leads to an implicit expression for Tc in the form: 
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where above, if I 
Eq.(g} reduces to: 
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electron bands are centered at 



-25 and 



-15, 
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where satisfies: Bf -z n B? + (tf-Jf)Bi+z n Jf = 0. We 
tested Eqs. (@J), JSJ in several cases: (1) at t?=t\i = and 
J 2 =0 we recovered the one-band model results reported 
in Ref. 0; (2) at < 2 =*i2=0, J 2 =0 and J x -» oo, the 
results of Ref. || are reproduced. Details on the calcu- 
lations above will be published elsewhere. From Eq. Q), 
the Tc contained in the Matsubara frequencies was ex- 
tracted numerically. Eq. J^J) predicts an increase in Tc 
only for the values of /x within the energy overlap interval. 

In Fig. Ufa), the DOS;(w) at J x /ti=25 and J 2 /t 2 =15 
is shown for several t U ' . Our Ji/ti are much larger than 
the value (J//^) m in ~ 1-4 that corresponds to the elec- 
tron and hole bands formation, and than the value J/t=8 
considered in the one-band model large enough to cap- 
ture properly the FM features of CMRs 0- Hence, the 



respectively (the hole bands, centered at w^=25 and 
cj^=15 due to the electron- hole symmetry, are not shown 
here for simplicity). For clarity, we plotted DOS2(w) with 
a reversed sign. Since | Ji/ii — Jz/tz] ^> 1, the bands oc- 
cupy different energy intervals. As a consequence, each 
band gives its own contribution to Tc (Fig. ^b)). At 
t U ' =0 the bands are fully decoupled, and the values of 
Tc match the results of the one-band model for all p's. 
However, if t U ' ^ 0, then the carriers are allowed to hop 
between the bands and, thus, they can belong simulta- 
neously to both bands. This leads to an increase in the 
effective number of interacting electrons of each active 
band. As shown in Fig. ^a), due to the transfer of elec- 
trons among bands, in DOS;(cj) a new region occupied 
by the interacting electrons builds up within the inter- 
val of energies occupied by the DOS;' (u). At t u '=ti, the 
effective number of interacting electrons within each en- 
ergy interval becomes twice as large. In consequence, the 
Tc for all p corresponding to (i within each energy inter- 
val, is almost twice larger than in the one-band case (Fig. 
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FIG. 1: (a) DMFT zero-temperature interacting DOS; for 
Ji=25 and J2=15 at different values of £12; (b) Tc vs. p 
for the cases shown in (a). Since the Hund's couplings are 
very different, the bands and critical temperatures retain the 
semicircular form; (c) Zero-temperature interacting DOS; at 
ii2=0.5, for the couplings indicated; (d) Tc vs. p for the 
cases shown in (c). Even if the overlap of bands is narrow, 
the exchange of electrons is strong enough to induce a robust 
Tc at p — 1 (black curve). In all frames t\=ti=\. 



^b)). By contrary, when the Ji's are similar the coupling 
of bands is strong enough to lead to a deviation from the 
semicircular form, especially when tw ~ \ Ji — J ; < | (Fig. 
H^c)). If the bands partially overlap, a hump develops in 
Tc at all values of p corresponding to \i within the energy 
overlap interval (Fig.^d)). 

In FigEIa), the total interacting DOS(u>) 0] at h=l, 
£2=1/3 ~ 0.33 is shown, for different t u > in the case 
when the bands fully overlap ( J\/J2=l)- The electron 
exchange effect does not only increase the effective num- 
ber of interacting electrons in each band, but also extends 
the energy region occupied by the bands. However, the 
total number of interacting electrons does not change, 
the bands being fully filled for p=2. At t 12 =0 (Eq. ©) 
the increase in Tc is due to the overlap of bands (green 
dashed curve in Fig. 2(b)). However, when t u > 7^ (see 
Eq. (g|) the exchange effect further boosts Tc (Fig-^Jb)). 

Monte Carlo results: The Hamiltonian Q was also 
studied numerically using the MC methods widely ap- 
plied to Mn-oxides in the limit J — > 00 0. Hence, 
the e g -spins are perfectly aligned with the i 2 g spin. The 
technique involves finding the eigenvalues of the Hamil- 
tonian matrix at each MC step corresponding to a newly 
updated set of localized spins. Although this substan- 
tially limits the size of the clusters being simulated, the 
results for small lattices are numerically exact and they 
allow for a direct comparison with DMFT. We have sim- 
ulated lattices of sizes 8x8 in 2D and 4x4x4 in 3D. The 
core spins are treated classically, while the treatment of 



FIG. 2: (a) DMFT zero-temperature total interacting DOS 
for different values of t\2- (b) Tc vs. p for the parameters 
indicated in (a). The red curve corresponds to the realistic 
set of hoppings for LaMnOa. In all frames ii=l, Ji = J2=15. 
The black-dot curve corresponds to the one-band case. 



the fermionic sector is exact. In all simulations 2 x 10 4 
MC steps were used, the first 10 4 been discarded in order 
to account properly for the thcrmalization of the random 
starting configuration. In finite dimensions, the hopping 
carries a direction index "a", that is tf v . In 2D (3D), 

3*f, = 1, tf,=V3tf„=x/3if, , = 



tf^-Vstf^-Vst 



i',r 



2,t v v v = 1, (tf d =tf r ^tf, l= Q, tf, ;/=4/3), in the x-, y- (and 

z-) directions, respectively The indexes I and I 

stand for the two active, x 2 — y 2 and 3z 2 — r 2 , orbitals. 
ii = l sets the energy unit. To find Tc we investigate the 
long-range spin-spin correlations: 
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where N is the total number of sites. Tc is the temper- 
ature for which S — > upon heating, at the maximum 
distance x max in the clusters considered. 

In Fig. [31 the results in 2D and 3D are displayed side- 
by-side for comparison. In panel (a), the spin-spin cor- 
relations at x=x max is shown for an 8x8 lattice, corre- 
sponding to different electron densities p=N e /2N, where 
A^e is the total number of electrons. The same in (d), but 
on a 4 3 lattice. Moreover, we investigated Tc vs. p at 
ti=l) t2 = l/3, an d different ii2- The results are in Fig. 
|3Jb) for 2D and in Fig. EJe) for 3D. The Tc is maximum 
at p=l and vanishes in the limits p^0 and p^2, with an 
overall semicircular form in the phase diagram, as seen 
in the frames (c) and (f). The one-band phase diagram 
(denoted lb) is also shown. Hence, as seen in (c) and 
(f), tu has a substantial effect in rising the To, which is 
in qualitative, and even in quantitative [l2T |. agreement 
with DMFT. As ii2 increases from to 1, the increase in 



4 



3-Gp=U 

3-ap=1/e 
p=1/2 

s-ep=5/4' 
7/4 
2 ■ 




FIG. 3: (a) MC spin-spin correlations at the maximum dis- 
tance 4\/2 vs. temperature for different values of p, at 
ti=l(= t), t2=l/3, and the t\2 indicated, using a 2D lat- 
tice of size 8x8. Errors are less than the symbol size, (b) 
Same as in (a), but changing ti2\ (c) Tc vs. p at different ti2 
for the same lattice used in (a). The black line corresponds 
to the one-band case. The red line shows the Tc obtained 
from (a), (d-f) Same as in (a-c), but using a 3D lattice of size 
4x4x4. In all frames Ji = J2 — > 00. 



T c can be as high as 100% in 2D and 60% in 3D. 

The finite-size effects are checked using different 
boundary conditions as well as simulating clusters of up 
to 12x12 in 2D and 5x5x5 in 3D. The results are in 
Fig. where the spin-spin correlations at x — x max (a) 
and the magnetization \M\ vs. T curves (b) are shown. 
The size effects are small, showing that our MC method 
detects the Tc accurately. 

Conclusion: We carried out the first study of a multi- 
band DE model applied to CMR using a powerful com- 
bination of DMFT and MC techniques. When two ac- 
tive bands are considered, the Tc is maximized at p—1. 
DMFT shows that the intcrband hopping leads to an in- 
crease in Tc at all p's, even if the electron bands do not 
occupy the same energy interval. This is due to the elec- 
tron exchange between the bands, which increases the ef- 
fective number of interacting electrons within each band. 
Both DMFT and MC indicate that, if the bands fully 
overlap, besides the increase of Tc due to the energy 
overlap, a further boost occurs when the intcrband hop- 
ping is turned on. The ideas developed here can be used 
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FIG. 4: MC finite-size effects studied in 2D dimensions, (a) 
Spin-spin correlations at the maximum allowed distance vs. 
temperature for the lattice sizes indicated, (b) Magnetization 
\M\ vs. temperature corresponding to the same parameters 
in (a). The non-zero value of \M\ at high temperatures is 
the asymptotic value 1/viV for a system of size N. 3D re- 
sults using 5x5x5 lattices are similar (not shown). While the 
spin-spin correlations do not show appreciable size effects, the 
magnetization results clearly do, their errors bars being max- 
imal within the critical region. 



to search for multiorbital FM materials with even higher 
Tc than currently known, once the interband hopping is 
tuned up. Our study can be extended to Diluted Mag- 
netic Semiconductors, with similar results expected |13j . 
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